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Abstract 

In the small dispersion limit, solutions to the Korteweg-de Vries equation de- 
velop an interval of fast oscillations after a certain time. We obtain a universal 
asymptotic expansion for the Korteweg-de Vries solution near the leading edge of 
the oscillatory zone up to second order corrections. This expansion involves the 
Hastings-McLeod solution of the Painleve II equation. We prove our results using 
the Riemann-Hilbert approach. 

1 Introduction and statement of result 

In recent works, critical behavior of solutions to Hamiltonian perturbations of hyper- 
bolic and elliptic systems has been studied. In particular the universal nature of this 
critical behavior is remarkable. For hyperbolic one- and two-component systems as 
well as for elliptic two-component systems, it has been conjectured by Dubrovin et 
al. [m [TU [16] that solutions can be described in a certain critical regime by special 
solutions to the Painleve I equation and its hierarchy. For several hyperbolic systems, 
two other critical regimes have been observed. One of those regimes, which is related 
to the second Painleve equation, is the subject of the present paper. 

As a prototype-example of a Hamiltonian perturbation of a hyperbolic one-com- 
ponent PDE, we consider the small dispersion limit where e ^ for the Korteweg-de 
Vries (KdV) equation 

ut + 6uux + e^Uxxx = 0, xGR, tGM+, e > 0. (1.1) 

We are interested in the KdV solution u{x, t, e) at time t > when starting at t = 
with initial data u{x, 0, e) = uo{x) which we assume to be negative, smooth, sufficiently 
decaying at ±00, and with a single negative hump. For a much wider class of initial 
data than the one we consider, it is known that the KdV solution exists at all positive 
times t. Up to a certain time the KdV solution u{x, t, e) can be approximated by the 
solution to the Cauchy problem of the dispersionless equation 

ut + (}uux = 0, u{x,0) = uo{x), (1-2) 

which is the Hopf equation. It is a classical fact that the solution to the Hopf equation 
reaches a point of gradient catastrophe at a finite time tc > 0, where the derivative of 
the solution blows up. After this time, the Hopf solution does not remain well-defined 
and can only be continued as a multi-valued function. The dispersive term e'^Uxxx in the 
KdV equation (jl.ip regularizes the gradient catastrophe for any e > 0. As a drawback, 
the KdV solution develops oscillations with wavelength 0{e) for t > tc, see Figured] 




Figure 1: Solution of the KdV equation with e = 10 ^, initial data uo{x) = — sech^(a;), 
and t = 0.4. 



Those oscillations were already described in the classical works of Lax and Lev- 
ermore [32j and investigated in more detail by Venakides [36], Tian [35j, and Deift, 
Venakides, and Zhou [HI US]. Under mild assumptions on the initial data, u{x, t, e) can 
be asymptotically described using Whitham equations [37J and elliptic ^-functions. Up 
to a fixed time after the time of gradient catastrophe, say for tc < t < T, there is one 
oscillatory interval [x~{t),x~^{t)]. For x £ {x~ (t) , x~^ (t)) (and x fixed), we then have 
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Here I3i{x,t), P2{x,t), ^^{x^t) evolve according to the Whitham equations, and 
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(1.4) 



K{sy K{s)' " /3i-/?3 

where K{s) and E{s) are the complete elliptic integrals of the first and second kind, 
K'{s) = K{VT^), and ^{z; t) is the Jacobi elliptic theta function given by the 
Fourier series 



For t > tc but X bounded away from the oscillatory interval [x {t),x^{t)], the asymp- 
totic behavior of the KdV solution is still determined by the Hopf equation. 

In the (x, t)-plane up to a certain time T > tc, we distinguish a 'Hopf region' and 
an oscillatory region which is cusp-shaped. The caustics separating the two regions 
consist of the leading edge, the trailing edge, and the point of gradient catastrophe 
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Figure 2: Cusp-shaped oscillatory region in the (x, t)-plane for t > tc- 

itself, see Figure [TJ Although we do not use this fact, we note that those borders are 
the curves on which the so-called Lax-Levermore maximizer [321 [36j is singular and the 
associated Riemann surface changes from genus to genus 1. The transitional asymp- 
totic description of the KdV solution u{x, t, e) for (x, t) approaching the caustics, has 
escaped investigation for a long time. Near the point and time of gradient catastro- 
phe, it was proved recently that the KdV solution can be approximated using a higher 
order Painleve I transcendent [5]. This supports the conjecture of Dubrovin [H], who 
claimed that this Painleve I-type behavior should be to a large extent universal, see 
also |15j . The behavior near the trailing edge has not been clearly understood so far, 
but we intend to come back to this problem in future work. Near the leading edge, 
multiple scale analysis and numerical results by Grava and Klein [28] indicate that the 
Hastings-McLeod solution to the Painleve II equation takes over the role of the higher 
order Painleve I transcendent. It is the aim of this paper to prove this rigorously for a 
large class of initial data. It is expected that the Hastings-McLeod solution to Painleve 
II arises universally near the leading edge for a whole class of perturbations to hyper- 
bolic equations, including (among others) the Camassa-Holm equation [H [29] and the 
de-focusing nonlinear Schrodinger equation. 

It is remarkable that the KdV solutions can be approximated by Painleve equations 
in the critical regimes described above. Those Painleve equations have appeared in 
many branches of pure and applied mathematics during the last decades, for an overview 
we refer to [22] . An application of particular interest for us is the appearance of the first 
and second Painleve equation in the theory of random matrices, see e.g. [2 H [T8 | [3| 16117] . 
In random matrix ensembles with certain unitary invariant probability measures, the 
eigenvalues accumulate on a finite union of intervals when the dimension of the matrices 
grows [8l [9l [To] . Including a parametric dependence to the probability measure for the 
matrices, it is possible to observe transitional regimes where the number of intervals in 
the spectrum changes. This can happen in three different ways. The first one occurs 
when two intervals merge to a single interval, and this situation shows remarkable 
similarities with the leading edge for the KdV equation. The second possibility, which 
should be related to the trailing edge for KdV, is that an interval in the spectrum 
shrinks and disappears afterwards. Finally it could happen that the merging of two 
intervals takes place simultaneously with the disappearance of one of the intervals. This 
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last phenomenon is related to the gradient catastrophe for KdV. 

The Painleve II equation is the following second order ODE in the complex domain, 

q"{s) = sq + 2q\s). (1.5) 

The special solution in which we are interested, is the Hastings-McLeod solution 
which is uniquely determined by the boundary conditions 



q{s) = V-s/2(l + o(l)), ass^-oo, (1.6) 

g(s) = Ai(s)(l + 0(1)), as s ^ +00, (1.7) 

where Ai(s) is the Airy function. Although any Painleve II solution has an infinite 
number of poles in the complex plane, it is known [31] that the Hastings-McLeod 
solution q{s) is smooth for all real values of s. 

1.1 Statement of result 

As already mentioned, we consider negative smooth initial data with only one local 
minimum and with sufficient decay at ±00. To be more precise, we impose the following 
conditions. 

Assumptions 1.1 (a) uo{x) is real analytic and has an analytic continuation to the 
complex plane in a domain of the form 

S = {z £C: \lmz\ < tan6'o|Rez|} U {z G C : |Imz| < a}, 

for some < Oq < tt/I and o" > 0, 

(b) uq{x) decays as \x\ -^00 in S such that 

uoix) = O (^j-^y d>0, xeS, (1.8) 

(c) U(){x) < and uq has a single local minimum at the point x = 0, with 

u'oiO) = 0, <(0) > 0, 
and Uq is normalized such that uq{0) = —1, 

Condition (b) is needed in order to apply the scattering transform [2] for the KdV 
equation. This transformation associates a reflection coefficient to the initial data uq. 
Condition (a) yields certain analyticity properties of the reflection coefficient. The 
presence of only one local minimum as imposed in condition (c) is convenient in order 
to have the simplest possible situation in the study of the semiclassical asymptotics for 
the reflection coefficient, while the position of the minimum and the minimal value of 
Uq can be chosen to be and —1, respectively, without loss of generality. 

The gradient catastrophe for the Hopf equation takes place at time 
max^giR;[-6^i'o(4)J 

Let us denote fL for the inverse function of the decreasing part of the initial data uq 
(the part to the left of 0), so that 

t, = -i^max^^[/i(0] = -IfM. (1.10) 
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The gradient catastrophe is generic if f'l'{uc) 7^ 0, and we will assume that this is the 
case for our initial data uq . It is known [35\ [30] that there exists a time T > tc such 
that for tc < t < T, the leading edge x~(t) is determined uniquely by the system of 
equations 

x-{t) = 6tuit) + fL{u{t)), (1.11) 
6t + e{v{t);uit)) =0, (1.12) 
dj{v{t);u{t)) =0, (1.13) 

with u{t) > v{t) and with 

Furthermore x~{t), u{t), and v{t) are smooth functions oft. Throughout the rest of the 
paper, whenever we refer to u, we mean by this the solution of the system p.ll|) - p.l3p 
for a given time tc < t < T, while we denote the solution of the KdV equation as 
u{x,t,e) and the solution of the Hopf equation as u{x,t). The equations (jl.lip - (jl.l3p 
correspond to the Whitham equations in the confluent case where 

/9i(x,t) = f32{x,t) = v{t), Psit) = u{t), 

see also (|1.3p . Note that the point of gradient catastrophe corresponds to /?i(x,t) = 
P2{x,t) = P3{x,t) (and thus u{t) = v{t)) and the trailing edge to the case where 
I3i{x,t) < (32{x,t) = (33{x,t). 

We will prove the following result. 

Theorem 1.2 Let uo{x) be initial data for the Cauchy problem of the KdV equa- 
tion satisfying the conditions given in Assumptions and furthermore such that 
the generic condition f'l'{uc) 7^ holds with Uc given by U.10\) . Then there exists a 
time T > tc such that the following is true for any fixed t such that tc < t < T . Take a 
double scaling limit where we let e ^ and at the same time we let x — > (t) in such 
a way that 

lim ''~^J^^ =X gM. (1.15) 

In this double scaling limit, the solution u{x,t,e) of the KdV equation with initial 
data Uq has the following asymptotic expansion, 

4e^/3 f@(x,t) i/s^ , 
u{x, t,e) = u -pr-q [s{x, t, e)\ cos h e ' fc)i(x, t, e) 

X - x~ 4e^/^ , , ,,0 9/9(x,t)\ ^, ^ , 

6t + c^l^(u — v) V £ / 

Here x~ and v < u (each of them depending on t) solve the system and the 

phase Q{x,t) is given by 

e(x,t) = 2V^7^(x-x-) + 2 / (/£(^) + 6t)7e^de (1.17) 



Furthermore 
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with 9 defined by 1^1. 14\ ), and q is the Hasting s-McLeod solution to the Painleve II 
equation. The correction to the phase Qi{x,t,e) takes the form 



Qi{x,t,e) 



.1/3 



q "^J 6d^29{v;u) A{u - v) 



^ s{x,t,ef / dlie{v]u) _ 3 ^ 2cy/V^ 



(1.19) 



4 \?,dl9{v-u) 2{u-v) m + f'^iu) 
where we used the notations 

q = Q{s), q' = q'{s), P = P{s) = -g'^is) - sq'^{s) + q'{sf, 
with s = s{x, t, e). 

Remark 1.3 Note that the leading order term in the expansion (jl.l6p of u{x,t,e) is 
given by u{t). The second term in (|1.16p is of order e^/^, while the further terms are 
of order e^/^. From the 0(e^/^)-term, we observe that u{x,t,e) develops oscillations of 
wavelength 0(e) at the leading edge. Zooming in at the leading edge, the amplitude of 
the oscillations is of order e^/^ and proportional to the Hastings-McLeod solution q. If 
we formally take the limit X — > — oo (so that x lies to the left of the leading edge), the 
terms with the oscillations disappear due to the exponential decay of q, see ()1.7p . We 
are then left with only two terms in (11.16p . which are the first two terms in the Taylor 
series of the Hopf solution u{x,t) near x^ . This is not surprising since we know that 
u{x, t, e) is approximately equal to u{x, t) outside the oscillatory interval. 

Remark 1.4 The expansion (11.16P is in agreement with the numerical results obtained 
in [28]. There is also a remarkable similarity between (jl.l6p and an asymptotic expan- 
sion obtained in [H [7j for the recurrence coefficients of orthogonal polynomials with 
respect to critical exponential weights. Those polynomials describe eigenvalue statis- 
tics in unitary random matrix ensembles. This can be interpreted as one instance of 
the more general analogy between Hamiltonian perturbations of hyperbolic equations 
and random matrices, which we discussed before. 

Remark 1.5 The time T in Theorem 11.21 should be chosen in such a way that the 
system (Il.lip - (jl.l3p is solvable for tc < t < T, and furthermore in such a way that 
Proposition 13.11 below holds true for tc < t < T. We do not go into further detail 
here, but we note for the interested reader that these conditions are equivalent to the 
regularity of the Lax-Levermore maximizer for tc < t < T and x = x~(t). 

Remark 1.6 It is unclear to what extent our result remains valid if the gradient catas- 
trophe is not generic, i.e. /^'(iic) = 0. For special choices of non-generic initial data, it 
is likely that the KdV asymptotics are described by higher members of the Painleve II 
hierarchy. 

The proof of our theorem relies on a Riemann-Hilbert (RH) problem which charac- 
terizes solutions to the KdV equation. This RH problem can be constructed using the 
direct scattering transform, while solving the RH problem corresponds to the inverse 
scattering transform. First of all in Section O we construct the RH problem and we col- 
lect some known asymptotic results about the reflection coefficient which is associated 
to the initial data uq via the direct scattering transform. Using those asymptotics, we 
can start performing the Deift/Zhou steepest descent method [13\ [TT\ [T2] in order to 
find asymptotics for the relevant RH problem. This is the content of Section [3l where 
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we will apply a number of transformations to the RH problem in order to have suit- 
able jump conditions. After a contour deformation, we will obtain a RH problem with 
jump matrices which decay uniformly to constant matrices when e — > 0, except near 
the special points u and v. We will then construct an outside parametrix away from u 
and V, and two local parametrices, one near u and one near v. The local parametrix 
near u will be constructed using the Airy function, and the one near v using special 
functions related to the Painleve II equation. After matching the three parametrices, 
we obtain uniform asymptotics for the solution to our RH problem in Section HI and 
consequently also for the KdV solution u{x^ t, e) when x approaches the leading edge 
of the oscillatory zone at an appropriate rate. 



2 RH problem for KdV and asymptotics for the reflection 
coefficient 

In this section we state the RH problem for the KdV equation, and we give a number 
of asymptotic results concerning the reflection coefficient for the Schrodinger equation 
with potential uq. The RH problem has been constructed in [Ml [2], and will be the 
starting point of our asymptotic analysis later on. The asymptotic results for the 
reflection coefficient were obtained in [33j and are crucial in order to perform a rigorous 
Deift/Zhou steepest descent analysis of the RH problem. A somewhat more detailed 
construction of the RH problem and the reflection coefficient can also be found in [5]. 

2.1 Construction of the RH problem 

The KdV equation (II. ip can be written in the Lax form using the operators 

o d , < , A 1 ^ ( f d d , ^ 

L = e — ^ + u[x, t,e), A = Ae — r + 3 u[x, t,e)— + —u{x, t, 
ax^ ax'^ \ ax ax 

The KdV equation can then indeed be written as 

L = [L,A], (2.1) 

where L = is the operator of multiplication by ut{x, t, e). A remarkable consequence 
of this Lax form is that the spectrum of L does not depend on time if u{x, t, e) evolves 



according to the KdV equation [25]. For negative initial data uq, it is known that the 
Schrodinger operator L with potential uq has no point spectrum 



Another important consequence of the Lax form (j2.ip for the KdV equation is 
that the solution u{x, t, e) can (for any x, t, e) be recovered from a RH boundary value 
problem, where one searches for a complex 2x2 matrix-valued function M which is 
analytic in C \ M and for which the boundary values on M are connected in a prescribed 
way. Suppose that M satisfies conditions of the following form. 

RH problem for M 

(a) M{X;x,t,e) is analytic for A G C\M. 

(b) M has continuous boundary values M+(A) and M_(A) when approaching A S 
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. \ {0} from above and below, and 

'0 1 
1 



M+(A) = M_(A)ai, ai=(° i), forA>0, 



with a given by 

a(A; X, t) = 4t(-A)2/2 + x(-A)^/^ (2.2) 

defined with the branches of (—A)" which are analytic in C\ [0, +oo) and positive 
for A < 0. Furthermore for fixed x, t, e, we impose M{X) to be bounded near 0. 

[ ' M 

(c) M{X;x,t,e) = {I + 0{X-^)) { as A ^ oo. 

For a suitable function r appearing in the jump matrices, the above RH problem is 
(uniquely) solvable, and the solution to the KdV equation with initial data uq can be 
recovered from the RH problem by the formula 

uix, t, e) = -2ie5^Mi,ii(a;, t, e), (2.3) 

where Mii{X;x,t,e) = 1 H ^' ^^ ^ ' — + 0(1/A) as A ^ oo. This is indeed known 

[2llllj to be the case if r(A; e) is the reflection coefficient from the left for the Schrodinger 
equation with potential uq, 

(P 

Lof = Xf, Lo = e^-^ + uq{x), 

see |19[ [2] for the definition and basic results about the reflection coefficient. The 
RH solution M can be expressed explicitly in terms of fundamental solutions to the 
eigenvalue equation Lf = Xf. As L has no point spectrum, those fundamental solutions 
are not L^-functions. 

2.2 Asymptotics for the reflection coefficient 

The reflection coefficient r(A; e) is independent of x (and t, since it is the reflection coef- 
ficient at time 0) but depends on the initial data uq via the direct scattering transform. 
Solving the Riemann-Hilbert problem corresponds to the inverse scattering transform. 
In general, there is no simple expression for r(A;e) in terms of uq. What is known, is 
that r(A; e) is analytic in a sector of the form vr — 29q < argA < vr (for some 6q > 0) 
if Assumptions 11.11 are satisfied [24j • Furthermore a number of asymptotic results for 
r(A; e) are known as e ^ 0. For A < —1, it is known that the reflection coefficient 
decays exponentially fast as e — > 0. On the interval (—1,0) on the other hand, r(A;e) 
oscillates rapidly for small e. To be more precise, the following holds as e ^ 0, 

r(A;e) =ze-TPW(l + 0(e)), for AG (-1,0), (2.4) 

r(A;e) = 0(6"?), for A G (-oo, -1), (2.5) 

where p(A) is given by 

/o(A) = /l(A) + / WMx) - A - V^]dx, (2.6) 
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and where denotes, as before, the inverse of the decreasing part of the initial data 

We need some more precise estimates for the reflection coefficient, in particular near 
— 1, where the transition takes place from oscillations to exponential decay. Because of 
the assumptions we imposed on the initial data uq, we can write 

/l(A) = -<7l(A)7aTT, (2.7) 

where ql is real analytic in a neighborhood of the interval [—1,0) and positive at — 1, 
and where \/A + 1 is analytic in C \ (— oo, —1] and positive for A > —1. We then have 

pW = -5l(A)V-A(A + 1)+ / WMx) - A - V^]dx. (2.8) 



Now it is easy to see that p is real analytic on (—1, 0) and that it has branch points at 
— 1 and 0. Observe that p extends analytically to the region 

1^+ := {A G C : -1 - 25 < Re A < -25, < Im A < 26} (2.9) 

for sufficiently small 6. Using the analyticity of the reflection coefficient in one can 
thus define an analytic function k{X) as follows, 

k{X; e) = -ir{X; e)eTpW, for A G n+. (2.10) 

In |33j . Ramond proved that, for any sufficiently small 6 > 0, we have the following as 
e-^0, 

K(A;e) = 1 + 0(6), uniformly for A E and Re A > -1 + 2(5, (2.11) 

r(A;e) = ©(e-f), for A < -1 - 25, c> 0. (2.12) 
In addition 

1 - |r(A;e)p ~ e-^^(^\ for -1 + 25 < A < 0, (2.13) 

where 

rfRW , 

t(A) = / ^/X-uo{x)dx, (2.14) 
J MX) 

and fji is the inverse of the increasing part of the initial data uq (the part to the right 
of the origin). Near —1, say for A G 1^+ and Re A < —1 + 25 (again for 5 sufficiently 
small), another result in [33] states that 

/^(A;e) = . ' (l + 0(g)), ase^O, (2.15) 
with given by 



vre 



where T is the Euler F-function. Since r maps the region {A G il+,ReA < —1 + 25} into 
the upper half of the complex plane, does not meet any poles of r(l/2 + z), and 

it follows that k(A; e) is bounded for A G $7+, Re A < —1 + 25. Moreover, if we delete a 
small semi-disk of radius 5/2 centered at —1 from it follows from the relation 

N{z) = 1 + 0{z^^), as z ^ CO, \ argzj < vr, 

together with (f2m) that 

k(A; e) = 1 + (D{e), uniformly for A G |A + 1| > f , as e ^ 0. (2.17) 
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3 Asymptotic analysis of the RH problem near the lead- 
ing edge 

From now on, we fix tc < ^ < for sufficiently small T > tc and, when convenient, we 
suppress t-dependence in our notations. We assume also that |x — x~| < 60 foi some 
sufficiently small > 0. Following the Deift/Zhou steepest descent method, we will 
apply a number of transformations to the RH problem for M. Here we follow similar 
lines as in [111 I12j . where a similar analysis was carried out away from the caustics 
in the (x,t)-plane, and as in [5], where a Riemann-Hilbert analysis was used near the 
point of gradient catastrophe. The most crucial new feature is the construction of local 
parametrices near the points u and v. 

3.1 Construction of the G-function 

Define G as follows 

G[\-x) = / —d'q, (3.1 

71" Ju \/V - u{v - A) 

where the choice of u is unimportant for now, as long as u is independent of x and 
— 1 < n = u{t) < 0. Later on, it will be important that we choose u according to 
(frTT]) - (frT3]) . In dni]), p is as in I^B), and a is given as before in UO^ . 



q(A; x) = 4t(-A)3/2 + x(-A)^/2_ (32) 

The square root in front of the integral in (13. ip is chosen to be analytic in C \ [u, +00) 
and positive for X < u, so that G is analytic in C \ [u, +00) as well. G has the following 
asymptotic behavior, 

G{X;x) = Gi{x)-^— + 0{X-'), as A ^ 00, (3.3) 
(-A)V^ 

where Gi is given by 

Gi(x) = i r^M^^d,. (3.4) 



Note that p{r]) is independent of x, so that 

^,G^{x) = -- r .f^dr^ = ^. (3.5) 



Ju \ V-u 2 

Using Cauchy's theorem and keeping track of the branch cuts for the square roots, it 
is easily verified that G has the properties 

G+(A) + G_(A)-2p(A) + 2a(A) =0, asAG(u,0), (3.6) 

G+(A) + G_(A) = 0, asA>0. (3.7) 

3.2 Transformation M ^ T 

We define 

T(A; X, e) = M(A; x, e)e-^^(^'^)'^^ (3.8) 

with (T3 = . T solves a RH problem with modified jump matrices compared 

to the ones for M. Using (|3.8p and (|3.6p - (|3.7p . we have the following RH problem for 
T. 
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RH problem for T 

(a) T is analytic in C \ M, 

(b) T+(A) = r_(A)wr(A), as A G M, with 



vt{X) = < 



1 



r(A)e" 



KG(A)+a{A))^ 



^_^(A)e-f(^W+-W) l-|r(A)|2 ^ 

^g-f(G+(A)-G-(A)) ^^^^^^p(A) 

-f(A)e-f (1 - |r(A)|2)e7(^+W-^-W) 



(c) r(A) ~ ( . . ^, — r- 1 as A — > oo. 

\i\/—X —i\/—X, 



as A > 0, 

as A G (— oo, u), 

as A G (u, 0), 



(3.9) 



as A — > oo. 



(3.10) 



(3.11) 



Using (I23D and ([331) -([SS]), we obtain 

u{x, t,e) = u — 2iedxTi^ii{x, t, e), 

where 

Tn(A;x,t,e) = 1 + + ©(A'^), 

Let us now define (/)(A; x) by 

(f){X; x) = G{X; x) — p{X) + a(A; x). 

This function is real for A G (— 1, -u) and it has branch points at —1 and u. Furthermore 
(p extends analytically to the region given by (|2.9p . By (|3.6p we have 

20+(A;x) = G+(A;x) - G_(A;a;), forAG('u,0). (3.12) 

The definition of (f) is helpful because it enables us to write the jump matrices for T in 
a more convenient form. Let us write (/'*(A) = (f){X) and k*(A) = k(A), with k given by 
()2.10p . Using (|3.9p together with (j3.12p and reality of p on (— 1, 0), we obtain 



vt{X) 



'^g-f 0+(A;x) 



i«;(A; e) 

m*(A;e) (1- |r(A)|2)eT<^+(^;^) 



as A G {u, 0). 



(3.13) 



Furthermore we also have 

1 



vriX) 



^iK*(A;e)e-T<^*(^;^) 



iK(A;e)eT<^(^;-)' 
l-|r(A)p ^ 



as A G [-1 - 5,u), (3.14) 



where for A G [—1 — (5,-1], 0(A), (/>*(A), k(A), and k*(A) should be interpreted as the 
boundary values where Im A \ (for the functions without star) and Im A /" (for the 
functions with star). The formula ()3.14p for A G (— 1,-u) follows directly from ()3.9p and 
the reality of (p. In order to obtain this formula as well for [—1 — (5,-1], note that the 
12-entry of vt (given by ()3.9p ) is analytic in while the 21-entry can be extended 
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to an analytic function in The same holds true for the right hand side of (|3.14|) . 
and by the identity theorem, we obtain (j3.14p for X G [—1 — 5, —1]. 

If condition (jl.lip holds, one verifies that (see [5, Lemma 3.2]) 

cPiX- x) = V^^{x -x-)+ J^fdO + Qt)y^l^d^. (3.15) 
In addition we have 

(/.'(A; x) = (x - x") - Vu - X[6t + e{X; u)], (3.16) 

2vn — A 

where 6{X; u) has been defined in (jl.l4p . If we let u = u{t) and v = v{t) be defined by 
(fTTT^ - ffOBI) . it follows that 

(j)'{v;x~) = 0, (f)"{v;x~) =0, (3.17) 



(P"'(v; x" ) = -V^^r^e"{v; u)=:c>0. (3.18) 
Proposition 3.1 There exists a time T > t^ such that for t^ < t < T , we have 
(j)'{X;x-) > 0, as X £ [-l,u)\{v}, 

Im0+(A;x-) < 0, asXe{u,0], (3.19) 

-t{X) + i((>+{X;x~) < 0, asXe{u,0]. 

Proof. For t = tc, it follows from (jl.lOp that the function /L(A)+6ic has a double zero 
at A = n = Uc (see also [5]) and that it is strictly negative on [—1,0] \ {uc}- It follows 
from (|1.14p and u{tc) = v{tc) = Uc that the same holds for the function 6tc + 6{X; u(tc))- 
For t slightly bigger than tc, because of the smoothness of u{t), 6t + 9{X; u{t)) can thus 
have at most two zeroes on [—1,0]. By ()3.17p . we know that both of those zeroes 
must coincide at the point v. Consequently we have that Qt + 9{X;u{t)) < for A S 
[—1,0) \ {v} for sufficiently small times after the time of gradient catastrophe. This 
implies the first inequality. We also have lm(j)'_^_{X;x~) < for X > u, and together 
with <j){u;x~) = 0, this implies the second inequality. 

For the last inequality, if r is defined as in (12.14p , it is straightforward to verify that 



3 f" /"^ 
-TiX) + iM^;x-) = -itiX-u)-2+ j ^./x^^f'^{S)di- j ^,/x^U'^{^)di, 

where fa is the inverse function of the increasing part of the initial data uq{x). The 
inequality follows directly because each of the three terms on the right hand side are 
negative. □ 

Remark 3.2 The smallest time T for which the above inequalities do not hold any 
longer, is the first time after the gradient catastrophe at which the Lax-Levermore 
maximizer becomes singular for x = x~{t). Our RH analysis is only valid at times for 
which (|3l9ll holds. 

3.3 Opening of the lens T ^ S 

Note that, by (|3.13p . we have the following factorization of the jump matrix for T, 
"^(^'"'^)= U.*(A;e)e-T^;(^-) 10 1 ' ^ ^ M " 
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Figure 3: The jump contour after opening of the lens 



(3.20) 



Because of this factorization, we can deform the jump contour to a lens-shaped contour 
as shown in Figure [331 Note that we start opening lenses at a point —1 — 5, and not at 
— 1 itself. This is a modification compared to the analysis in [T2] and is necessary 
to have jump matrices for S which are uniformly close to constant matrices. We will 
determine the precise shape of the lenses later on. 

As (/)+ can be extended analytically to 0+ and to 0+, we can define S as follows, 
\ -iK(A;e)eT<^(^;^)\ 



S{X) 



r(A) 



r(A)( 







1 



1 







IK 



*(A;e)e-T'^*(^;^) 1 



inside the lenses in the upper half plane, 

inside the lenses in the lower half plane, 
elsewhere. 



Now the RH problem for S takes the following form. 

RH problem for S 

(a) S" is analytic in C \ S5, 

(b) 5+(A) = S-{\)vs for A G S5, with 



1 iK{\)e 

^iK*(A;e)e-T<^*(A;a=) ij ^ 

iR{X;e) (1- |r(A;e)|2)eT'/'+(^;^) 



on Si 



on S2 



(c) Six) 



vriX), 



1 



as A G [u, 0), 
as A G (-00, -1- 6)U (0, +00). 



1 



as A ^ 00. 



J\/—X —i\/^X, 
Since T(A) = S{X) for large A, we have, using (|3.1Up . 

u{x, t,e) = u — 2iedxSi^ii{x, t, e), 



(3.21) 



(3.22) 
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where 

SniX;x,t,e) ^ 1 ^ -^i.iiM. e) +0(A-^), as A ^ oo. (3.23) 

In the foUowing proposition, we show that the jump matrices for S tend uniformly to 
constant matrices as e — > 0, except near the special points u and v. 

Proposition 3.3 For any fixed neighborhoods U of u and V of v, there exists Sq > 
such that for \x — x~\ < 5q, the jump matrices vs{X) decay uniformly to the matrix 
(A; X, e) on S5 \ (Z^ U V U {0}) ife^O, with 



«(°°)(A) 



0"!, on (0, +00); 

icji, on (u, 0), (3.24) 



/, elsewhere. 



Proof. For A > 0, our claim is trivial. For A < —1 — 5, we have that vs = vt and 
using ()3.9p , exponential decay to the identity matrix simply follows from (I2.12p together 
with reality of G and a. For A G (it, 0) \ if x = x~, we have uniform exponential 

decay of the 11-entry of vs (see (I3.2ip ) using the second inequality in Proposition 13.11 
and of the 22-entry using the third inequality in Proposition 13.11 together with (j2.13p . 
Furthermore if x is close to x~ , the first term in ()3.15p can only spoil the inequality in 
U and V. The off-diagonal entries tend uniformly fast to i using (j2.1ip . 

Using the first equality in Proposition 13. H it follows from the Cauchy-Riemann 
conditions that, again for x close to x~ , lm(p{X;x~) > on the upper parts of the lens 
if Re A > —1 + (5, and lm(p*{X; x~) < on the lower parts of the lens if Re A > —1 + 5, 
at least this is true if we choose the lenses sufficiently close to the real line. Again the 
inequalities remain valid for x close to x~ except in U and V. 

Near —1, we have by (j3.15p that 

Im(f>{\;x) = {x-x-)lmVu^+lm j {fHO + Qt)y^C - Xd^. 

It is then straightforward to check that (for x sufficiently close to x~) Im(j)(X;x) < 
on the upper parts of the lens near —1, and Im0*(A;x) > on the lower parts of 
the lens near —1, including the point —1 — 5 itself. Together with ()2.17p . this implies 
convergence of vs on the lenses near —1 as well. □ 

3.4 Outside parametrix 

Ignoring neighborhoods U and V of n and v, and ignoring uniformly small jumps, we 
are left with the following RH problem. 

RH problem for 

(a) : C \ [n, +00) C2x2 is analytic, 

(b) p(°°) satisfies the following jump conditions on (u,+oo), 

pjoo) ^ p^oo)^^^ (0,+oo), (3.25) 

pM^^^M^^^ on(7x,0), (3.26) 
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(c) has the following behavior as A ^ oo, 

This RH problem is solved by 




(3.28) 



and we have the following behavior at infinity, 

p(-)(A) = (/+^^3 + 0(A-2))(^.^_l^^/2 -i(-A)V2)' asA-oo. (3.29) 

The outside parametrix will determine the leading order asymptotic behavior of u{x, t, e) 
in the double scaling limit. However, in order to do a rigorous asymptotic analysis and 
to obtain correction terms in the asymptotic expansion for u{x,t,€) (which is our main 
goal), we need to construct local parametrices in neighborhoods of the special points u 
and V. 



3.5 Local parametrix near u using the Airy function 

We will construct a parametrix Pu, in a disk U surrounding u, which has approximately 
(as e ^ 0) the same jumps as S and which matches with the outside parametrix 
in the double scaling limit where e — > and simultaneously x — > x~ in such a way that 
X — x~ = ©(e^/^). To be precise, the aim of this section is to construct a parametrix 
Pu = -fu(A; X, e) satisfying the following conditions. 



RH problem for Pu 

(a) Pu is analytic mU \ E5. 

(b) Pu satisfies the jump conditions 



^-0{A;x) 



P„,+(A) = P„,_(A) J ), asAGEin^/, (3.30) 



P„,+(A) = Pu,- (A) ( .g^0.(A;x) ^ ) ' as A E S2 n Z^, (3.31) 



Pu,+{X) = Pu,-W[ . Q 1 , as A > u. (3.32) 

(c) As e — > and x ^ x~ in such a way that x — x~ = 0{e^^^), Pu matches with 
in the following way, 

P„(A) = p(°°)(A) (/ + 0(1)) , for A G dU. (3.33) 
3.5.1 Model Airy RH problem 

We will use a model RH problem using the Airy function which is not precisely the 
same but very similar to the one used in [8j and in many other works. It involves 
functions yo,yi,y2 solving the Airy differential equation and given by 

= 2/j(C) = e— Ai(e— C), J = 0,1, 2, 
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where Ai is the Any function. Let 



^i(C) = -zV2vre-x- { "^^ eX-, (3.34) 



^2(0 = -i\/27re-x-3 ( ] gx-a^ (3.35) 

"2/2 2/1^ 



^3(0 



,V2^e-x-3('yo (3.36) 
2/1/ 



Each of the functions Ai, A2, and are analytic in the entire complex plane. Using 
the well-known identity yo + 2/i + 2/2 = 0, it is easy to verify that 

^1(0 = ^2(0 (J ly (3.37) 
^2(0 = ^3(0 (J J), (3.38) 

Ai(C) = ^3(C)(j 0). (3.39) 

Furthermore the asymptotic expansion at infinity for the Airy function shows that 

A.(C) = e-x-3(_^)-^^ ^j + o e^(-C)^^^-3^ (3.40) 

uniformly as — > 00 in the sector Sj, j = 1, 2, 3, with 

5,={CGC::^7r<argC<^vr}, j = 1,2,3. (3.41) 

The fractional powers in (|3.4U|) are principal branches, so (— C)^ is analytic for < 
aigC < 27r and positive for arg^ = vr, and is given by ()3.28p . 

3.5.2 Construction of the parametrix 

We will now construct explicitly in terms of A. Define, in the neighborhood U of u, 
C(A) such that 

:^i-C{X)f' = ^{X;x-). (3.42) 

It then follows from p.lSp that is a conformal mapping in a sufficiently small neigh- 
borhood U of u, with 

C'in) = {-f'du)-6tf''>0. (3.43) 
Next define a function s(A; x) by 

s(A; j;)(-C(A))^/2 ^ 0(A;x) - 0(A;x") = {x - x-)Vi^^, (3.44) 
so that s is analytic in U with 

s{u;x) = {-fL{u)-6ty^^^{x-x'). (3.45) 
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Now we can define the parametrix P as follows, 

P„(A;x,e) =E(A;e)Aj(e-2/3c(A)+e-2/3s(A;2;))e-^<^(^;^)-3, (3.46) 

where j = 1 in the region for A outside the lens in the upper half plane, j = 2 inside the 
lens-shaped region, and j = 3 outside the lens in the lower half plane. The pre-factor 
E given by 

E{X; e) = p(°°)(A)A^-n-e"^/^C(A))^"^eX-3. (3.47) 

It is easily checked using (|3.28p that E is analytic in U. 

Using the jump conditions (|3.37|) - (|3.39p for A, it follows from (|3.46|) that Pu indeed 
satisfies the required jump conditions. More subtle is the matching of Pu with 
Let e ^ and at the same time x ^ x~ in such a way that x — x~ = 0(e~^/^). Then 
also s(A; x) = ©(e^/^), and we have by (102]) and i^M), 

-| (-^-^/^C(A) - e-^/hW^" = + for A e dU. 

Inserting the asymptotics (|3.4U|) for A into (j3.46p now yields the following matching of 
Pu with 



P4A)=P(°°)(A) 



is{X;xy s{X;x) , 



for A G dU, (3.48) 



in the double scaling limit. This matching will contribute to the asymptotics for the 
KdV solution u{x, t, e) later on. 

3.6 Local parametrix near v using ^-functions for Painleve II 

We will now construct the local parametrix in a sufficiently small disk V surrounding 
V. This construction shows similarities with the construction of the local parametrices 
in [HIT]. The parametrix P„ = P„(A;x,e) should satisfy the following conditions. 

RH problem for P„ 

(a) Pv is analytic in V \ S^. 

(b) Pu satisfies the jump conditions 



Pu,+ {X) = Pu,-W \ , asAESinV, (3.49) 

Pu,+ {X) = P,,_(A) (^.^^r(A;x) fj^ as A e E2 n V. (3.50) 

(c) In the double scaling limit where e ^ and simultaneously a; — > x~ in such a 
way that x — x~ = 0{e'^^^), Pu satisfies the matching 

P4A) =p(°°)(A)(/ + o(l)), forAeaV. (3.51) 
3.6.1 Model RH problem related to Painleve II 

We will construct the parametrix explicitly using a model RH problem for the so-called 
^'-functions associated to the Painleve II equation, and in particular related to the 
Hastings-McLeod solution of Painleve II. This RH problem appeared several times in 
the literature, see e.g. [22l [20l [Tt [3| [6] . 
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RH problem for ^ 

(a) ^' : C \ (Fi U C2><2 is analytic, where 

Ti = e*iR+ Ue^^M+, Fa = e~*iM+ U e'^^M" 
and FijFa are both oriented towards the right. 

(b) ^ satisfies the following jump conditions on Fi and F2, 



*+(C;s) = *_(C;s)(} J), asCGFi, (3.52) 

^+(C;s) = ^-(C;s)fn asCeF2. (3.53) 



.0 1 ^ 

(c) ^'(C; s)e^(t^'+"^)'^3 = J + 0(C"^) as C ^ 00. 

This RH problem is related to the Painleve II equation in the following way. Con- 
sider the linear system 

The compatibility condition ^'^^ = ^'^^ of this system implies that q solves the Painleve 
II equation (jl.5p and that y{s) = q'{s), and therefore this system is called the Lax pair 
associated to Painleve II. If we let q be the Hastings-McLeod solution to Painleve II, 
characterized by (ll.6p - (ll.7p . and y = q' , the RH problem for ^ can be solved using 
fundamental solutions to the Lax pair, see [201 122j . The RH problem is solvable for 
every complex value of s at which q{s) does not have a pole. It is known [31] that 
q is free of poles on the real line, so that the RH problem for ^ is solvable for s in 
a neighborhood of M. Condition (c) of the RH problem holds moreover uniformly for 
s in compact subsets of C \ where V is the set of poles for q, and the asymptotic 
condition can be strengthened as follows, 

.)e^(t«^+^^)- = I + ^ + + 0{C% as C - 00, (3.54) 



with 



where 

p{s) = -q\s) - q\s)s + q'isf, b{s) = ^ + lq{s)pis), (3.56) 

and a(s) is unimportant for us. 

In order to create a RH problem with jumps which model the jumps which are 
necessary for the parametrix, we insert an additional parameter uj and define 

We then have the following RH problem. 
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RH problem for ^ 

(a) $ is analytic in C \ (Fi U r2 U M). 

(b) $ satisfies the following jump conditions, 

1 



cI>+(C) = 1>-(C) ( ,^2^.,-2.(|c^+<) J > on r^. (3.59) 



(c) As C — > cx), we have 



$(C; = e-'^''-^ i^I + + + 0(C^') J e*""^ (3.60) 

where 

^'^'^ - 2l \-iq{s) p{s) ) ' ^'^'^ - 2i \h{s) a{s) ) ^^'^^^ 

This behavior is valid uniformly in cu, and uniformly for s in compact subsets of 
C\V. 

3.6.2 Construction of the parametrix 

We construct explicitly in terms of the function Let us take of the following 
form, 

P,{X;x,e) = P^^\X)'^ (^e-'/'CW;e-'/h{X;x),-e-'^{v;x)) , (3.62) 

where s{X;x) is analytic in V, and where CW is a confer mal mapping in V which maps 
V to 0. We can then specify the precise shape of the lenses S5 near v by requiring that 
C(5]) C Fi U so that P has its jumps on S5. 

We will now determine ("(A) and s(A; x) in such a way that 

cPiX; x) - cl^iv; x) = ^C(A)' + s(A; x)C(A). (3.63) 

If this condition is satisfied and if we let 

uj{x,€) = — (j){v;x), (3.64) 

the jump conditions (I3.49p - (l3.50p for P^ are valid using the jump conditions for <I> and 
analyticity of near v. We can organize and s such that (|3.63p is valid, if we first 
define ( by 

cl){X;x~)-^{v;x~) = ^CW\ (3-65) 
so that using (j3.17p . we have 



19 



where c has been defined in (j3.18p and 9{v; u) in (|1.14|) . Next define s(A; x) by 
s{X;x)({X) = (j){X; x)—(j){X; x~ )+(j){v; x~ )—(j){v; x) = {x—x^){V u — X—\/u — v), (3.67) 
so that we have 

X — x~ X — x~ 



s{v; x) 



s'{v; x) = s{v] x) 



2('{v)\/u — V c^/^\/u — V ' 
1 lC"{v) 



4(m -v) 2 C'{v) 



Let us now take the double scahng limit where e ^ and at the same time x 
such a way that e^^/^(x — x~) ^ X G M. This means that 



-.X. 



(3.68) 
(3.69) 

■ x~ in 

(3.70) 



C^/^y/u — V 

In this double scaling limit where in addition we let A approach v, we also have that 

lime-2/=^s(A;x) =s + 0{X-v). 

This implies that e~'^/'^s(A; x) does not meet any poles of q and lies in compact subsets 
of C \ "P for A in a sufficiently small neighborhood V of v. Summing up (j3.65p and 
()3.67p now gives indeed (j3.63p . 

If A is at a fixed small distance away of v, for example on 9V, we have that 
e^^/^("(A) — > oo if e ^ 0. Now we can use the asymptotic expansion ()3.60p for 
Inserting this into ()3.62p . we obtain the following matching in the double scaling limit 
for A E dV, 



P(A;x,e) =pM(A)e-^'^'^3 

1/3 2/3 

/ + ±^$i(e~2/3^(^. + l_^^^,-2/z,^x- x)) + 0{e) 



3.7 Final transformation S ^ R 



(3.71) 



Define 



R{X; x, e) 



f5(A)P„(A)-i, 
5(A)P.(A)-i, 
[5(A)P(-)(A)-1, 



in U, 
in V, 

elsewhere. 



(3.72) 



Using the jump conditions for S and the ones for the parametrices, one verifies that 
the jump matrices for R are uniformly close to the identity matrix in the double scaling 
limit. This will imply that R itself is also uniformly close to the identity matrix. In 
particular we have RH conditions for R as follows. 



RH problem for R 

(a) R is analytic in C \ S/j, where S/j = S5 U dU U dV. Here we choose the clockwise 
orientation for dU and dV. 
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(b) R satisfies the jump condition R+[z) = R-{z)vr{z) on S/j, where vr has the 
fohowing asymptotics in the double scahng hmit, 



vr{z) 



' I + Vi{X]x,e)€^/^ + V2{\]x,e)€^/^ + 0{e), as A e dVUdU, 
l + 0{e), as A G (ZYU V) nEfi, 

/ + 0(6-^1/^), asAGSR\Z7UV. 

(3.73) 



(c) R{X) = 1 + 0{X-') as A ^ oo. 

The uniform decay of the jump matrix outside U DV rehes on Proposition 13.31 and the 
jumps for the outside parametrix. The jumps inside U and V are easily verified to be 
I + 0{e) using (j2.17p . the jumps for S, and the jumps for P. The most important 
jumps (i.e. the only ones which will contribute to the expansion (11.161) for u{x, t, e)) are 
the ones on dU and dV. Using (|3.48p and (j3.28p one verifies that, for A G dU, we have 



(3.74) 



4e4/3(-C(A))i/2 V(w-A)V2 

and for A G dV we have using (|3.7ip . 

Vi{X-x,e) = p(°^)(A)e-^^"^$i(e-2/3s(^.a.))ei^-3p(oo)(^)-i 
_ p{e-^/^s{X-x)) ( (?/-A)-i/2^ 



2C(A) \-{u-\Y'^ 

q{e~'^/^s{X;x)) ( cos2w(a;,e) (u - A)-!/^ sin 2w(x, e) 



2C(A) A)i/2sin2w(x,e) -cos2w(a;,e) ^"^'^^^ 

Note that s{X;x) = ©(e^/^), so that Vi remains bounded in the double scaling limit. 
Let us now take a look at V2. On dU, we find the following by (I3.48P and (|3.28p . 

and on 9V, 

V2{X-x,e) = p(°°)(A)M2P(°°)(A)-i 

_ a[e-^l^s[X;x)) 
2<(A)2 

h{e-'^/^s{X-x)) ( -sin2a;(x,e) {u - co^2uj(x,t)\ ^^^^^ 



2C(A)2 A)i/2cos2w(x,e) sin2w(a;,e) 

Note that V\ and V2 can be extended to meromorphic functions in Z// U V, where V\ has 
simple poles at u, f, and V2 has a simple pole at u and one of order two at i;. Y\ and 
V2 remain bounded in the double scaling limit, but they do depend on e. 
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4 Asymptotics for R and the KdV solution u{x, t, e) 

Because the jump matrix for R has an asymptotic expansion in powers of e^^^, it follows 
[SI [To] that R itself has a similar expansion (uniform in A) in the double scaling limit, 

R{X; x,e) = I + e^/^R^'^'^ (A; x, e) + e^/^R^^'^ (A; x, e) + 0{e). (4.1) 

We will compute large A asymptotics for R^^^ and R^'^^ using the following proposition. 

Proposition 4.1 As A oo, we have 

i?W(A;x,e) = i(Res(yi;t;) +Res(yi;M)) + 0(A-2), (4.2) 
A 

R^^\X;x,e) = \ iRes{W;v) + Res{W;u)) + 0{X-^), (4.3) 
A 

where 

W{X) = /?W(A)yi(A) + V2{X). (4.4) 
Proof. Combining (|3.73p with (j4.ip gives us 

+ gi/3^(i) (A) + ( A)) ^ = + eV3^(i) ^ ^2/ 

X (^/ + yi(A)ei/2 + V2{X)e^^'^^ + 0(e), for A G 5V U OZ^. 

Collecting the terms of order 0{e^^^) and of order 0{e^^^) leads to the identities 

i?W(A) = i?L^)(A) + Vi{X), for A G 5^/ U 9V, 

i?(2)(A) = R^^\X) + i?i^)(A)Vi(A) + V2(A), for XedUU dV. 

Together with analyticity of R^^^ and R^'^\ and the asymptotic condition 

i?(^)(A) ^0, as A ^ oo, j = 1, 2, 

the above jump conditions constitute a uniquely solvable RH problem for R^^^ and R^'^\ 
One verifies that, since Vi has simple poles at u,v, 

rW is given by 

/A^Re«(^i;^) + A^Res(^i;^) forAGC\(zYuv), 

\^Res(yi;?;) + ^Res(yi;n)-yi(A) for A G U V. 
Since V2 has a simple pole at u and a double pole at v, we have 

) f^A) forA.C\(«UV). 
\f(A)-VF(A) forAGZ^UV, 

with 

F(A) = -^Res(VK; v) + -^Res(T^; u) + ^ ,^ Res((A - v)W; v). 
X — V X — u (A — t>)^ 

This imphes (g^]) and (gS]). □ 
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As A — > oo, R can be written as follows, 

R{X;x,e) = I+^^^^ + 0{\-^). (4.7) 
A 

Using the fact that S{X) = R{X)P^°°\X) for A away from u and v, using (|3.22|) and 
(|3.29p we get the following identity, 

u{x, t,e) = u — 2edxRi,i2{x, t, e). (4.8) 

As a consequence of (I4.2|) - (l4.3p . it follows that 

Ri = {Res{Vi;v) +Res{Vi;u)) e^^^ + {Res{W;v) +Res{W;u)) e^^^ + 0{e) (4.9) 

in the double scaling limit. For notational convenience, let us write q := q{e~'^^^s{v; x)), 
and similarly for a,b,p,q' in the remaining part of the paper. In what follows, we 
carefully analyze the residues in (j4.9p and the x-derivatives of their 12-entries. 



Proposition 4.2 In the double scaling limit, we have 

X — x~ 



e^/3a,Res(yi,i2; 



2(6t + /[(n))' 



(\-.\ ^4/3 n /T. . „N _ 2gcos(2a;) q^ 0/3 g^sin(2u;) 0/3 

2^u — vc^l-^\%t + Ji[u)) 
/■ N 5/3 /Tx. X g2cos(4w)e2/3 2gsin(2u;)e2/3 /26 



C2 1 



6 4(n — i;) 



u — V 



3L r Q 
C2 3 



12 8{u-v) 



{x - x-)h~^/^ ^ c^/^x - x-)2e~4/3 



c2/3 (u-v) 4 Vu^(6t + (n) ) y 

""^ e{v;u) 



with UJ given by 6'^[ ), c 6y ([ITT^), a'^fi C2 = 4S 

Proof. (i) The first equation follows directly from (|3.74p using (j3.43p and (j3.45p . 

(ii) Note first, using (^3M^ and H^Mh that 

dxuj{x; e) = -\/u - ve"\ dxs{v; x) = — . ] 

c'-i'^yu — V 

Using these identities and taking the derivative of the residue of the 12-entry in 
(|3.75p . we find equation (ii) using ()3.66p . 

(iii) By we have 

Res(T^;y)i2 = i?S(y)Res(Fi; y)i2 + 42^(y)Res(Fi; y)22 + Res(F2; y)i2, (4.10) 

both ioT y = u and y = v. For y = u, we have that Res(Vi; n)22 and Res(V2; y)i2 
vanish. The only nonzero term is thus the first one. Since we know that r\^^ (n) = 
Res(Vi;t;)ii/(u - v) by (g^]), we conclude from i!^7I^ and ([37751) that 

q s'm(2uj)s(u: x)'^ oi-i , , , . 

d.ResiW; u)u = \ / =^^ ' + 0{e). (4.11) 

4Vn-<(v)VC'(^i) 

Inserting the formulas (j3.43p . (j3.45p . and (j3.66p . we obtain (iii). 



23 



(iv) Taking a look at (j4.1Up for y = v, all three terms are nonzero. The third term at 
the right hand side of (|4.10p gives by (|3.77p . 



Res(y2;v] 



12 



cos(2a;) 



2C'{v)^^/^r^ \ Civ) ' 2{u-v) 

For the first two terms in (I4.10p . we use ()4.5p and obtain 



cos(2w) 



(4.12) 



pC"{v) 



2C{v)y/u^ \ 2C'{v) 

sin(2a;) /gC"(^^) 



p's'{v;x)e '^^^ 



P 



2{u - v) 



2C'{v)y/u^ \ 2C'{v) 



's'(t;;x)e-2/3 



2{u-v)_ 



4(u-i;)vrM 



Now we need to use the fact that 2b = q' + qp and that p' = —q^ (see (|3.56p 
and (jl.Sp ). After a straightforward but somewhat lengthy calculation, using the 
above equations together with (j3.74p - (j3.75p and ()4.10p . we find 



q^ sin(4a;) 



gcos(2u;) 

16C'(i')2(u-i')3/2 ' ^CWWu^ 



+ 



2bC"{v) P+k 
qC'{v) u — V 



s'{v; x)s{v; x)e 



-4/3 



C{v)s{u;xf 

2,Ar^y^ciu)' 



-4/3 



Taking derivatives and using (ITOD . (lOSD . (^M), and (IMH1)-(IM9|), (iv) is ob- 
tained directly. 

□ 

Now we can insert the above estimates into (j4.9p and (|4.8p . However we need to 
be careful with the derivative of the C'(e)-term in ()4.9p . On dU, one can check that 
this term is of order 0(e) by computing the next term in the expansion (j3.48p . On dV, 
using (j3.7ip . we have a contribution of order 0{uje) = 0{1). Also on the other parts of 
the contour, we have a contribution of order 0{e). Taking this into account, we obtain 



u{x, t, e) = u{t) - 4— ^jyj— e 



_ gcos(2cj) _ 4g^ sin'^(2u;) 2/3 , a: - 



c^/^{u — v) 



6t + f^iu) 



4gsin(2w)e2/3 / 2b 
\ 1 



.2/3 



C2 



6 A{u-v) 



£2 
12 



(u — v) 



c^/^iu — v) 



c^l^{x^x~ft-^l^ 
'2Vir^(6i + /[(u)) 



+ 0(e). (4.13) 



Observe that 2a; = - \Q{x,t), with a; defined by (f3:M]) and G by (fTTTD . After a brief 
calculation, we can now write (j4.13p in the form (|1.16p . so that Theorem 11.21 is proven. 
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